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In this study, the rapidity distribution, collective flows, and nuclear stopping power in 1°7Au +19" Au 
collisions at intermediate energies were investigated using the ultrarelativistic quantum molecular dynamics 
(UrQMD) model with GEMINI++ code. The UrQMD model was adopted to simulate the dynamic evolution of 
heavy-ion collisions, whereas the GEMINI++ code was used to simulate the decay of primary fragments pro- 
duced by UrQMD. The calculated results were compared with the INDRA and FOPI experimental data. It was 
found that the rapidity distribution, collective flows, and nuclear stopping power were affected to a certain extent 
by the decay of primary fragments, especially at lower beam energies. Furthermore, the experimental data of 
the collective flows and nuclear stopping power at the investigated beam energies were better reproduced when 


the sequential decay effect was included. 
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I. INTRODUCTION 


Understanding the properties of dense nuclear matter is 
a main topic of the basic research on heavy-ion collisions 
(HICs), as it is the only way of creating dense nuclear mat- 
ter in terrestrial laboratories [1-5]. However, the typical time 
scale for the appearance of the created dense nuclear mat- 
ter is extremely short (less than or of the order of a few 
fm/c) [6]; thus, direct measurement of its properties is cur- 
rently impossible. Usually, the fundamental properties of 
dense nuclear matter are inferred by comparing experimen- 
tal data and transport model simulations [7-11]. As pow- 
erful tools for investigating the dynamics of nonequilibrium 
systems at intermediate energies, two typical transport mod- 
els, the Boltzmann—Uehling—Uhlenbeck (BUU)-type [1] and 
quantum molecular dynamics (QMD)-type models [2] and 
their updated versions, have been extensively used for sim- 
ulating HICs to extract the structure of the initial nuclei and 
the properties of the created dense nuclear matter, such as the 
nuclear equation of state, the nuclear symmetry energy, and 
the in-medium nucleon-nucleon cross section [12-17]. 

Both the BUU- and QMD-type models mainly include 
three components: the initialization of the projectile and tar- 
get nuclei, mean-field potential for each particle, and two- 
body scattering (collision term). Different techniques are typ- 
ically used to treat the three components in different models. 
To establish a theoretical systematic error that quantifies the 
model dependence of transport predictions and disentangle 
the causes of different predictions, the transport model eval- 
uation project (TMEP) has been recently performed [18-22]. 
After simulating the dynamical process of HICs within the 
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transport models, an afterburner is usually chosen for the for- 
mation of fragments, which can then be used to construct ob- 
servables, and further compared with the corresponding ex- 
perimental data to extract interesting information regarding 
the nuclear matter. However, the formed fragments are typ- 
ically excited and cannot be experimentally detected. In the 
experiments, the detected particles underwent long-term de- 
excitation. Hence, greater focus should be placed on the ef- 
fect of de-excitation in the treatment of afterburners, espe- 
cially for studying HICs at Fermi energies where the beam 
energy is comparable to the exciting energy [23-25]. 


The underestimation of the yield of light fragments (e.g., 
3H, ?He and “He) in HICs at intermediate energies is a long- 
standing problem in transport model simulations, and numer- 
ous techniques have been developed and adopted to overcome 
it [26], including the statistical multifragmentation model 
(SMM) [27, 28], the statistical evaporation model (HIVAP) 
[29], the statistical model GEMINI [30], and the Simulated 
Annealing Clusterization Algorithm (SACA) [31, 32]. From 
our previous studies [33-35], adopting the conventional phase 
space coalescence model “Minimum Spanning Tree” (MST) 
[36-38]and disregarding the decay of the excited primary 
fragments from the Ultra-relativistic Quantum Molecular Dy- 
namics (UrQMD) model simulation, the experimental data of 
the collective flow of protons and deuterons at INDRA and 
FOPI energies can be reasonably reproduced, whereas those 
of a particles are different from the experimental data. To im- 
prove the corresponding collective flow distribution and ex- 
plore the influence of sequential decay on the observables in 
the HICs at intermediate energies, the statistical code GEM- 
INI++ was employed to describe the decay of primary frag- 
ments. 


In this study, the rapidity distributions and collective flows 
of free protons and light clusters, as well as the nuclear stop- 
ping power, were calculated based on the UrQMD model with 
and without considering sequential decay. The remainder of 


this paper is organized as follows. In Sect. II, the relevant de- 
scriptions of the transport model, statistical model, and the 
observables are provided. The calculation results are pre- 
sented and discussed in Sect. II. Finally, a summary and 
outlook are provided in Sect. IV. 


Il. MODEL DESCRIPTION AND OBSERVABLES 
A. Dynamical model: UrQMD 


The UrQMD model [39-42] is a typical transport model 
used for microscopic many-body nonequilibrium dynamics. 
In this model, each nucleon is represented by a Gaussian wave 
packet of a certain width L. Empirically, L = 2 fm? are chosen 
to simulate Au+Au collisions. The coordinates and momen- 
tum of each nucleon are propagated using Hamilton’s equa- 
tions of motion. The total Hamiltonian of the system (H) 
consists of the kinetic energy T and the effective interaction 
potential energy U, which includes the Skyrme potential en- 
ergy Up, Coulomb energy Ucou; and momentum dependent 
potential energy Uma: 


(H) =T + Up + Ucout + Uma- (1) 


To study the HICs at intermediate energies, the Skyrme en- 
ergy density functional was introduced in the same manner as 
that in the improved quantum molecular dynamics (ImQMD) 
model[43, 44]. The local and momentum-dependent potential 
energies can be written as Up ma = f up,madr where 
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Here, ô = (Pn — Pp) / (Pn + Pp) denotes the isospin asym- 
metry defined by the neutron(p,,) and proton (pp) densities. 
In the present work, the normal nuclear matter density po = 
0.16 fm~3, a = —393 MeV, 6 = 320 MeV, y = 1.14, gsur 
= 19.5 MeV fm?, gsuriso = —11.3 MeV fm?, Agym = 20.4 
MeV, Bsym = 10.8 MeV, Cyym = —9.3 MeV, n = 1.3, tma = 
1.57 MeV, and ama=500 (GeV/c)~? are chosen, correspond- 
ing to a soft and momentum-dependent (SM) equation of state 
with the incompressibility o = 200 MeV and the symmetry- 
energy slope parameter L = 80.95 MeV, which are in their 
commonly accepted regions [42, 45-49]. 

In this work, the in-medium nucleon-nucleon (NN) elas- 
tic cross section is considered as the result of a density- 
and momentum-dependent medium correction factor F(p, p) 
times the free nucleon-nucleon cross section o&°, which is 
available from the experimental data. 
= F(p, p) * ONN» (4) 
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(5) 
The parameters fo = 1, A = 1/6, ¢ = 1/3, po = 0.3 GeV/c and 
k = 8 were adopted, with pyy being the momentum in the 
two-nucleon center-of-mass frame. More detailed studies re- 
garding the effects of F(p, p) parameters on various observ- 
ables can be found in Refs. [41, 42, 50]. In addition to the in- 
medium N N cross section, Pauli blocking plays an important 
role in determining the rate of collisions. The Pauli blocking 
treatment is the same as that described in Refs. [19, 41]; it 
involves two steps. First, the phase space densities f; and fj 
of the two outgoing particles for each NN collision are cal- 
culated. 
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(6) 
Here, o? denotes the width parameter of the wave packet. o? 
= 2 fm? is used in this work. k represents nucleons with the 
same type around the outgoing nucleon 2 (or j). The follow- 
ing two criteria were simultaneously considered: 
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Prock = 1 — (1 — fj — fj) < £. (8) 


Here, rig and Pig denote the relative distance and momen- 
tum between nucleons 7 and k. The above conditions are also 
considered for nucleon j. The symbol € denotes a random 
number between 0 and 1. 

In addition, the isospin-dependent minimum spanning tree 
(iso-MST) algorithm was used to construct the clusters [36— 
38]. Nucleon pairs with relative distances smaller than Ro 
and relative momenta smaller than P = 0.25 GeV/c are con- 
sidered to be bound in a fragment. Here, RP = 2.8 fm and 
Ro” = Rj” = 3.8 fm are used for proton-proton and neutron- 
neutron (neutron-proton) pairs, respectively. A fairly good 
agreement between the recently published experimental data 
(including collective flows, nuclear stopping power, and Han- 
bury—Brown-Twiss interferometry) and UrQMD model cal- 
culations was achieved with appropriate choices of the above 
parameters [35, 42, 45, 51, 52]. 


B. Statistical model: GEMINI++ 


The statistical GEMINI++ decay model [30, 53] is a Monte 
Carlo method to simulate the decay of excited fragment, 
which includes light particle evaporation, symmetric and 
asymmetric fission, and all possible binary-decay modes. It 
has been widely used to treat the de-excitation of the primary 
fragments in studies of nuclear reactions at low and interme- 
diate energies [23, 25, 54-60] and high energies [61-63]. The 
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parameters of the shell-smoothed level density in GEMINI++ 
were set to their default values, which were ko = 7.3 MeV and 
Koo = 12 MeV. A more detailed discussion of the GEMINI++ 
parameters can be found in Refs. [53, 60, 64]. 

There are four inputs for the GEMINI++ code (with de- 
fault parameter settings): the mass number A, charge num- 
ber Z, excitation energy E*, and angular moment L of the 
primary fragment. After the UrQMD transport and iso-MST 
coalescence processes, the excitation energy of a fragment is 
calculated as follows: 

E* = Epa — Bima > (9) 
where | Ocoee is the binding energy of the primary fragment, 
that is, the sum of the potential and kinetic energies of all 
nucleons that belong to a fragment. E80" is the binding 
energy of the ground state obtained from nuclear data table 
AME2020 [65, 66]. When the excitation energy of the frag- 
ment was less than or equal to zero, it was regarded as a bound 
state and did not decay further. The angular momentum of the 
primary fragments was calculated using classical mechanics 
[56]. 


L= 7 x hi, (10) 
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where 7; and p; are the coordinate and momentum vectors, 
respectively, of the i-th nucleon in the primary fragments in 
the c.m. frame of the fragment. The total angular momen- 
tum is the sum of all the nucleons in the primary fragments. 
We verified that the results remained almost unchanged if the 
contribution of the angular momentum in GEMINI++ was not 
included. 


C. Observables 


In this study, we focus on collective flows and nuclear stop- 
ping power, which are commonly used observables, and pro- 
pose investigating the properties of dense nuclear matter. The 
directed (vı) and elliptic (v2) flows are defined as [67] 
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Here, pz, Py are the two components of the transverse mo- 


mentum p; = ,/pz + p3}. Angle brackets indicate the aver- 


age of all considered particles from all events. The nuclear 
stopping power Rpg, which characterizes the transparency of 
the colliding nuclei, can be defined as [68], 


2 
B= 3B (12) 


where E; (Ej) is the c.m. transverse (parallel) energy, 
and the sum runs over all considered particles. Another 
quantity vartl is also widely used to measure the stopping 
power, which was proposed by the FOPI Collaboration [69], 


is widely used to measure stopping power. It is defined as the 
ratio of the variances of the transverse to those of the longitu- 
dinal rapidity distribution and is expressed as 
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where (y2) and (y2) are the variances in the rapidity distri- 
butions of the particles in the x and z directions, respectively. 


vartl = 


MI. RESULTS AND DISCUSSIONS 
A. Rapidity distribution 


The rapidity distributions of the hydrogen (Z = 1) and he- 
lium isotopes (He and tHe) in 197 Au +19” Au collisions at 
Fiap = 40 MeV/nucleon and impact parameter b = 5.5-7.5 
fm are shown in Fig. 1. The left, middle and right panels 
are the calculated results when the switching time of UrQMD 
model was set to 200 fm/c, 250 fm/c and 300 fm/c, respec- 
tively. When considering the decay of the excited primary 
fragments, the yields of Z = 1 particles are slightly increased, 
whereas those of *He and *He clusters are clearly increased. 
The yields of He and tHe are approximately 1.5 and 6 times 
greater than those obtained from the simulations without con- 
sidering sequential decay. It is known that the yields of 
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Fig. 1. Rapidity distributions of Z = 1 [panels (al-a3)], ?He [pan- 
els (b1-b3)] and He [panels (cl-c3)] particles as functions of the 
reduced rapidities(yo = yz/Ypro) from peripheral (b = 5.5-7.5 fm) 
197 Ay +197 Au collisions at Ejay = 40 MeV/nucleon at different 
reaction times. The solid and dashed lines represent the results ob- 
tained from UrQMD simulations without and with considering the 
decay of the primary fragments by using GEMINI++, respectively 
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Fig. 2. Reduced rapidity yo distributions of the directed (elliptic) 
flow vı (v2) from semi-central (peripheral) °” Au +197 Au col- 
lisions at Fiab = 40 MeV/nucleon for Z = 1 particles. The blue 
circles denote the calculated results without primary fragments de- 
cay, whereas the olive down-triangles represent the calculated results 
with primary fragments de-excited 


light clusters such as 3He and *He obtained using QMD-type 
codes are significantly smaller than experimentally measured 
yields. As discussed in our previous work [44], the yield of 
3He calculated using the UrQMD model was approximately 
three times lower than that of the experimental data. The in- 
clusion of GEMINI++ improved the description of light clus- 
ter production to some extent; however, the yields of light 
clusters were still underestimated. Additionally, the inclusion 
of GEMINI++ may have led to an overestimation of the yield 
of free nucleons. To solve these problems further, additional 
issues should be considered, such as the spin degree of free- 
dom and the dynamical production process [7]. 


B. Collective flows 


The collective flows of Z = 1 particles as functions of yo 
in 197 Au +19” Au collisions at Ejay = 40 MeV/nucleon are 
shown in Fig. 2 (directed flow vı [left panel] and elliptic flow 
v2 [right panel]). The simulations are stopped at 200 fm/c. 
Clear differences are observed in the distributions of the col- 
lective flows between simulations with and without consid- 
ering the decay of the primary fragments, especially in the 
case of elliptic flow v2 from peripheral (b = 5.5-7.5 fm) col- 
lisions. From Fig. | the difference in the yields of Z = 1 par- 
ticles between simulations with and without considering se- 
quential decay is due to the particles produced from the decay 
of excited primary fragments, which have larger mass. These 
heavier excited primary fragments are usually produced at the 
target/projectile rapidity; thus, the directed flow at the tar- 
get/projectile rapidity is significantly affected by sequential 
decay. Increasing the impact parameter results in more large- 
mass primary fragments being produced and the effects of 
sequential decay on the collective flows becoming more evi- 
dent. 

To quantitatively evaluate the influence of sequential de- 
cay on the collective flows, as shown in Fig. 3 the vı slope 
(top panels) and v2 (bottom panels) at mid-rapidity for Z = 1 
particles (left panels) and free protons (right panels) are cal- 
culated at different stopping times and compared with the ex- 
perimental data. The slope of directed flow v1; and elliptic 


E,4p=90 MeV/nucleon — E,,,=400 MeV/nucleon 


rre rr 0.6 
0.1 | (a) b=2-5.5 (mee AMD ainil (2) D=3.4-6.0 fm FOPI| 
INDRA U?0.8 
= 0,0} tsa, rs > Pessa AS a 
> L | > 
-0.1> v--__) ---¥ -- I 10.4 
-0.2 fi fi fi 1 1 f fi A fi 1 
| (b) b=5.5-7.5 fm |_(d) b=3.4-6.0 fm +-0.04 
0.12}+ J 
Kaie) - J -0.06 g 
Sng. Shh A pg fb JL —+— SF > 
> 0.09 ULLAL. + -0.08 
0.06 C 1 1 | 1 1 fi Y2 Yorol<0.4 ] 
150 200 250 300150 200 250 300 
time (fm/c) time (fm/c) 
Fig. 3. Directed flow slope vi; and elliptic flow veo from 


197 Ay +19” Au collisions at Ea, = 50 MeV/nucleon (for hydro- 
gen isotopes (Z = 1)) and 400 MeV/nucleon (for free protons with 
uzo > 0.8) with different stopping times. The INDRA and FOPI ex- 
perimental data taken from Refs. [69-71] are indicated by shaded 
bands 


flow v9 at mid-rapidity are extracted by assuming v1 (yo) = 
Vio +011 Yo + V13 Ys and V2(Yo) = V20 + V22 ` YE +V24° YH 
in the range of |yo| < 0.4, in the same way as done in the ex- 
perimental report [69-71]. The left panels show the results 
from the simulations at Fiab = 50 MeV/nucleon, whereas the 
right panels show the results from the simulations at Fj,, 
= 400 MeV/nucleon. There is an obvious difference in v11 
and v29 between the simulations with and without consider- 
ing the decay of the primary fragments at F),, = 50 MeV/nu- 
cleon. By considering the decay of the primary fragments, 
the value of v;; decreases, while that of va9 increases, and 
both are closer to the experimental data. This can be under- 
stood from the fact that both directed and elliptic flows are 
stronger for large-mass fragments [33, 69, 70]. The decay 
of fragments is isotropic within the fragments c.m. in the 
statistical model. Thus, the light nuclei produced from the 
decay of primary fragments retain information about the pri- 
mary fragments. Furthermore, at Fiab = 50 MeV/nucleon, 
v11 first decreases with time and then saturates (increases) af- 
ter 200 fm/c in the case without (with) the sequential decay 
effect. This can be explained by the competition between se- 
quential decay and dynamical effects. At approximately t = 
150 fm/c, the nucleons are located in a low-density environ- 
ment, and the final state interactions are attractive; thus, a de- 
crease v11 is observed. When considering GEMINI++ at t = 
300 fm/c, there are only a few large-mass fragments that can 
emit hydrogen isotopes, and the effect of sequential decay on 
V11 gradually weakens. A comparison of the sequential de- 
cay effects at Fiap = 50 MeV/nucleon with those at Fiab = 
400 MeV/nucleon shows that the sequential decay effects are 
more pronounced at lower beam energies. This is understand- 
able because more large-mass fragments, which contribute to 
the production of light clusters, are formed at lower beam en- 
ergies. 

Fig. 4 shows the beam energy dependence of vj; (top 
panel) and vq (bottom panel) for semi-central 197 Au+1!97 Au 
collisions. The value of v11 from the simulations without con- 
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Fig. 4. Beam energy dependence of directed flow slope v11 (panel 
a) and elliptic flow v20 (panel b) from semi-central ° Au +19" Au 
collisions. As shown, the solid stars represent the data from INDRA 
Collaboration [70] for Z = 1 particles, and the open stars are the 
results from FOPI Collaboration [69, 71] for free protons with uto 
> 0.8 cut. The lines with different symbols are calculated with and 
without considering the sequential decay. Where error bars are not 
shown, they fall within the symbols 


sidering sequential decay is higher than that of the experimen- 
tal data at lower energies. When sequential decay is included, 
V11 is driven down and the vo is pulled up at low energies, 
which are consistent with the results shown in Fig. 3. Again, 
both vıı and v2ọ are considerably closer to the experiment 
data for all investigated energies. Moreover, the effects of se- 
quential decay on the collective flows decrease with increas- 
ing beam energy. With increasing beam energy, the reaction 
becomes more violent, and the multiplicities of large-mass 
primary fragments decrease. 

For a systematic investigation of the effects of the sequen- 
tial decay on the collective flows, the transverse 4-velocities 
uto dependence of the vı of light charged particles in semi- 
central !97Au +197 Au collisions at Ej, = 120, 150, 250 
MeV/nucleon are calculated with and without considering the 
decay of primary fragments, and compared with the experi- 
mental data [69], as shown in Fig. 5. The effects of sequential 
decay on vı (uzo) of particles with mass number A < 4 are 
relatively weak; it only affects vı at low uto. By consider- 
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Fig. 5. Directed flow vı of protons (squares), deuterons (circles), 
A = 3 clusters (triangles), and a-particles (rhombuses) versus uzo in 
197 Au+*°7 Au collisions at Fia» = 120, 150, 250 MeV/nucleon with 
b = 3.4-6.0 fm. Results from calculations with and without sequen- 
tial decay are represented by solid and open symbols, the FOPI ex- 
perimental data (red symbols) are taken from Ref. [69]. The shaded 
areas represent the corresponding error bars 


ing the sequential decay effect, vı of a particles is obviously 
influenced and can reproduce the experimental data reason- 
ably well at higher uo. This is because after the decay of 
excited primary fragments is considered, the multiplicity of a 
particles is enhanced, as shown in Fig. 1. The a particles pro- 
duced by the decay of the excited primary fragments have a 
large flow effect because they inherit the flow information of 
the excited primary fragments. The remaining discrepancies 
in uzo-dependent vı of light clusters may be due to simplifi- 
cations in the initial wave function of particles (nucleons and 
possible clusters) and quantum effects in two-body collisions, 
as well as the lack of dynamical production process, which 
requires further study [7]. 


C. Nuclear stopping power 


The nuclear stopping power characterizes the transparency 
of the colliding nuclei and can provide insight into the rate of 
equilibration of the colliding system. Thus, it is meaningful to 
explore the effects of sequential decay on the nuclear stopping 
power. The yield distributions of deuterons and *He clusters 
as functions of the reduced longitudinal (y,/yp,-o) and trans- 
verse (Yx/Ypro) rapidities for central 1°’ Au +19" Au colli- 
sions at Fiap = 40 MeV/nucleon with and without considering 
the sequential decay are shown in Fig. 6. The value of calcu- 
lated Rpg are shown in each panels, and they are slightly de- 
creased; that is, the nuclear stopping power is depressed when 
considering the sequential decay. Notably, the values of Rg 
are slightly influenced by sequential decay; however, the yield 
spectrum of each particle is obviously influenced, especially 
for He clusters. In conjunction with Fig. 1, it can be con- 
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Fig. 6. Yield distributions of 7H and “He particles as functions of the 
reduced longitudinal and transverse rapidities from central (b = 0-2 
fm) 197 Au +97 Au collisions at E,ap = 40 MeV/nucleon. Results 
from calculations with and without the sequential decay are shown 
in the left and right panels. The corresponding values for nuclear 
stopping power Rg are also indicated in each panel 


cluded that sequential decay has a strong effect on the parti- 
cle yield and yield distribution; however, the nuclear stopping 
power is only weakly influenced because the de-excitation of 
the primary fragments in the GEMINI++ code is isotropic. 
The nuclear stopping power Rg (Erab = 40 MeV/nucleon, 
left panel) and vartl (Ejay = 400 MeV/nucleon, right panel) 
for free protons from central (b = 0-2 fm) 19 Au +197 Au 
collisions calculated at different stopping times are shown in 
Fig. 7. Similar to the collective flows shown in Fig. 3, the 
effects of sequential decay on the nuclear stopping power at 
lower energies are more obvious than those at higher energies. 
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Fig. 7. Rg and vartl for free protons from central (b = 0-2 fm) 
197 Au +197 Au collisions with Fia» = 40 and 400 MeV/nucleon at 
different stopping times. The symbols sets are the same as Fig. 1. 
The shaded bands are the corresponding INDRA and FOPI experi- 
mental data taken from Refs. [72, 73]. 


T e 
1.0 F Au+Au free protons F 1.0 
b=0-2 fm t 
0.9 : E 0.9 
ot! 0.8 i T J Lose 
> 
Be [ 
0.7 } 0.7 
:  —O—UrQMD 7 
: =W -with GEMINI++ Í 
eo k INDRA! v FOPI pee 
40 100 400 


Elab (MeV/nucleon) 


Fig. 8. Beam-energy dependence of Re (INDRA) and vartl (FOPI) 
for free protons from central (b = 0-2 fm) 197 Au +197 Au collisions. 
The symbols sets are the same as Fig. 5. Calculations with and with- 
out sequential decay are compared with the INDRA and FOPI ex- 
perimental data, which are taken from Refs. [72, 73], respectively 


In all cases, the results were saturated and were close to or 
covered by the experimental data above 200 fm/c. 

Fig. 8 shows the degree of nuclear stopping (Rg or vartl) 
in the central 197Au +197 Au collisions as a function of the 
beam energy. It is clear that both Rg and vartl decrease 
and are much closer to the experimental data when the GEM- 
INI++ code is applied, and the difference in the values of the 
nuclear stopping power gradually disappears with increasing 
beam energy. With the same impact parameter (centrality), 
more violent reactions occurred at higher beam energies, and 
less heavily excited primary fragments were produced. The 
difference with free protons is that they are produced from 
the decay of heavy excited primary fragments, which usually 
have a weak stopping power (smaller values of Rp and vartl) 
[73, 74]. As a result, the nuclear stopping power was slightly 
reduced by including the GEMINI++ code. 


IV. SUMMARY AND OUTLOOK 


In summary, the effects of the sequential decay of ex- 
cited primary fragments on the rapidity distribution, collec- 
tive flows of light nuclei, and nuclear stopping power in 
Au+Au collisions at intermediate energies were investigated. 
Primary fragments were produced using the UrQMD model, 
and sequential decay of the excited primary fragments was 
performed using the GEMINI++ code. It was observed that 
the sequential decay of the excited primary fragments have an 
obvious influence on the rapidity distribution and collective 
flows but relatively weakly affect the stopping power. This is 
because of the memory effect; that is, the light particles pro- 
duced from the sequential decay of the excited primary frag- 
ments inherit the collective properties of the excited primary 
fragments. Furthermore, the sequential decay effects grad- 
ually decreased with increasing beam energy because fewer 
large-mass fragments that can emit light particles were pro- 
duced at higher beam energies. More importantly, the ability 


to reproduce the relevant experimental data at the investigated 
energies was improved to a certain extent by the inclusion of 
sequential decay effects. 


These results are meaningful for extracting the properties 
of nuclear matter by comparing transport model calculations 
with experimental data. The effects of secondary sequential 
decay on the observables, such as the N/Z ratios and elliptic 
flow ratios of neutrons vs. hydrogen isotopes, will be dis- 
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